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Abstract 

Thermally excited phonon spectra of 2D Yukawa solids and liquids in the presence of an external mag- 
netic field are studied using computer simulations. Special attention is paid to the variation of wave spectra 
in terms of several key parameters, such as the strength of coupling, the screening parameter, and the in- 
tensity of the magnetic field. In addition, comparisons are made with several analytical theories, including 
random-phase approximation, quasi-localized charge approximation, and harmonic approximation, and the 
validity of those theories is discussed in the present context. 
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I. INTRODUCTION 



There has been a great deal of interest in studying collective dynamics of two-dimensional 
(2D) strongly coupled Coulomb systems (SCCSs) in the presence of an external magnetic field 
II, 2. 3, 4, siy,!?!. Early studies were mainly concerned with the 2D electron gas and/or liquid 
systems, and were largely motivated by their potential applications in the semiconductor industry, 
e.g., for metal-oxide- semiconductor (MOS) structures llL S |3|, |4|] . Recently, interests have turned 
to a slightly different system, i.e., the so-called strongly coupled complex (dusty) plasmas (SCDPs) 

Complex (dusty) plasma is a suspension of micron-sized charged dust particles in a weakly 
ionized plasma with electrons, ions, and neutral atoms or molecules [|8|,|91]. Therein, dust particles 
acquire a few thousand electron charges by absorbing surrounding electrons and ions, and con- 
sequently interact with each other via a dynamically screened Coulomb potential lllOn while un- 
dergoing Brownian motion due primarily to frequent collisions with the neutral molecules. When 
the interaction potential energy between charged dust particles significantly exceeds their kinetic 
energy, they become strongly coupled and can form ordered structures comprising liquid and solid 
states. Since the motion of charged dust particles in complex (dusty) plasmas can be directly 
observed in real time by using a video camera, such systems have been generally regarded as a 
promising model system to study many phenomena occurring in solids, liquids and other SCCSs 
at the kinetic level. 

The 2D SCDP has become particularly favorable in recent laboratory experiments, probably 
because complications due to the ion-wake effect [Illl.ll2n can be avoided in 2D systems, and thus 
the interaction between dust particles can be well approximated by the Debye-Hiickel/Yukawa 
potential 11131 Il4l]. </)(r) = (Q^/r) exp (— r/A/?), where Q is the charge on each particle, r is the 
inter-particle distance, and A/j is the screening length. Such a system can be fully characterized 
in equilibrium by the screening parameter, k = a/ A/?, and the Coulomb coupling strength, T = 
Q'^/ {aksT), where ksT is the temperature of the system and a is the Wigner-Seitz radius, given 
by a = (vm)"^/^ with n being the areal density of particles. Of great interest in 2D dusty plasmas 
are their collective and dynamical properties, such as the longitudinal and transverse wave modes, 
which have been studied extensively over the past decade in experiment s 1115 . 



theories 1121 . 



22, 
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25 



27 



280, and numerical simulations 
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27 
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28h (see, in 



particular, Ref. [|27l] for a review of recent developments). 
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In particular, there has been a growing interest recently in studying the collective dynamics of 



SCDPs in an external magnetic field ||29 



31 



32i 



33 



3411 ■ However, most of those studies were 



theoretical. For example, Uchida et al. [|29ll derived a wave dispersion relation for a 2D Yukawa 
lattice in a perpendicular magnetic field based on the harmonic approximation (HA) and compared 
it with the results obtained through a Molecular Dynamics (MD) simulation in a perfect triangu- 
lar lattice. Jiang et al. OOll studied the propagation of Mach waves in a 2D magnetized Yukawa 
system based on the random-phase approximation (RPA), where the effect of strong coupling was 
neglected. That work was subsequently extended to include the effect of strong coupling Bill 



24, 



2711 ■ Moreover, Sun et al. 



by adopting the quasilocalized charge approximation (QLCA) [jSL 
[I32I1 examined the influence of an external magnetic field on a dust lattice wave in a plasma crystal. 
Very recently, Farokhi et al. [13 311 derived analytically a wave dispersion relation in a 2D hexagonal 
Yukawa lattice in an external magnetic field with arbitrary direction. 



Although many theoretical proposals have arisen from the above listed studies [129 



3311 and experimentlists have also attempted to study the dynamics of magnetized dusty plasmas 



[|34|l . exploring these proposals in the 



technical constraints [|29L 



32l 



33L 



aboratory does not seem to be feasible yet due to many 



34|l . In this situation, it is desirable to perform numerical 



experiments that can validate the existing theories and that can serve as a guide for future labora- 
tory experiments. Other than the above-mentioned MD simulation by Uchida et al ^il^ . we have 
found only two related MD simulations in open literature: one studying static structures of 2D 
magnetized dusty plasmas by Liu et a/.[|350, and the other studying collective modes in a Coulomb 
system (2D strongly correlated electrons) by Ranganathan and Johnson [IgI, 7]. Given that, to the 
best of our knowledge, there are no computer simulations studying collective modes in magnetized 
Yukawa systems that cover both the liquid and solid states, we perform the numerical experiments 
here to fill the gap. 

First, we obtain thermally excited phonon spectra from direct simulations of 2D Yukawa sys- 
tems in a perpendicular magnetic field. The resulting dispersion relations are then compared with 
those derived in the above mentioned analytical theories for different values of the screening pa- 
rameter and the Coulomb coupling strength. The rest of the paper is organized as follows. In Sec. 
II, we provide details of our simulations of 2D Yukawa solids and liquids with perpendicular mag- 
netic field. In Sec. Ill, analytical formulas are given for dispersion relations in the RPA, QLCA 
and HA, and compared with simulations. A discussion of the results is given in Sec. IV, followed 
by a brief conclusion in Sec. V. 
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II. SIMULATION 



A. Algorithm 



370, in 



Our simulation is performed by using the Brownian Dynamics (BD) method 
which we assume that numerous and frequent collisions with the neutral-gas molecules give 
rise to a Brownian motion of each dust particle in the plasma. These collisions give rise to 
the well-known Epstein drag, which is quantified by the friction coefficient, or damping rate 7. 
We therefore need to introduce one additional parameter, the damping coefficient 7/co'pd, where 
^pd = [2Q^/(ma^)]^^^ is the nominal dusty plasma frequency, enabling us to fully characterize 
the system 113 811 . 

In our work, = 4000 particles are simulated in a square with periodic boundary conditions. 
Charged dust particles interact with each other via pairwise Yukawa potentials, while undergoing 
Brownian motions. The simulation follows the trajectory of each dust particle by numerically 
integrating and solving the corresponding Langevin equation and its integral Soi 4ol] . 

4v = -7v + — F + — Qv X B + A{t), 
at m mc 

where, as usual, v and r are the velocity and position of a dust particle, respectively, and F is 
the systematic (deterministic) force coming from interactions with all other dust particles in the 
system. Furthermore, B = {0, 0, B} is the external magnetic field, c is the speed of light in vac- 
uum, and A(t) is Brownian acceleration, respectively. Note that the friction coefficient 7 and the 
Brownian acceleration A(t) represent complementing effects of the same sub-scale phenomenon: 
numerous collisions of a Brownian particle with the neutral-gas molecules in the system. While 
the former represents an average effect of these collisions, the latter represents fluctuations due to 
the discreteness of collisions and is generally assumed to be well represented by delta-correlated 
Gaussian white noise. The friction coefficient and the Brownian acceleration are both related to 
the medium temperature through a fluctuation-dissipation theorem. We adopt here a Gear-like 
Predictor-Corrector method for our BD simulations, which was used successfully in simulating 



the shock wave propagation [|41,|42|], the heat conduction [j43l], and diffusion processes in SCDPs 



[44]. Detailed algorithms and procedures for performing these integrations can be found in Ref. 



4211 . The advantages of using the BD method are as follows. First, it brings the simulation 
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closer to real dusty plasma experiments by taking into account both the Brownian motion and the 
damping effect in a self-consistent manner. Second, it simplifies the simulation as no external 
thermostat is needed. 

Initially, charged dust particles are randomly placed in the square. The system comes to an 
equilibrium after a period proportional to I/7, typically around IO/7, except under conditions 
with extremely strong coupling and large magnetic field, in which case relaxation of the system 
becomes extremely slow. One would expect a faster relaxation of the system (and therefore a faster 
simulation) at higher damping rates. However, we choose here 7 = 10^^ x cupii to minimize the 
damping effect on collective modes, thus allowing us to concentrate on the effects of the magnetic 
field. We note that the effect of neutral-gas damping on collective modes in unmagnetized SDCPs 
has been discussed in detail elsewhere [281. 



B. Calculation of the wave spectra 



Wave spectra are determined from the current-current correlation functions. We follow the 
definition given in Ref. [|45ll . in which the current density of the system is 



1 ^ 

j(r,t) = -^5Zv,(t)5[r-r,(t)], 

1=1 



(2) 



while its Fourier transform is given by 



1 ^ 



(3) 



i=l 



where rj(t) and Vj(t) are the position and velocity of the ith particle at time t, respectively, and 
the index a labels Cartesian components x and y. Assuming that waves propagate along the x 
direction, i.e., q = {q, 0}, one obtains the current-current correlation functions 

^i!(g,^) = (j:(g,i)jx(g,o)), 

J±{q,t) = {j;{q,t)jyiq,0)), 

for the directions parallel and perpendicular to the wave propagation direction, respectively. One 
recalls that, without the magnetic field, these functions correspond to the longitudinal and trans- 



verse modes [|28 



450, respectively. Here the superscript * denotes the complex conjugate while 



the angular bracket denotes the ensemble average. The wave spectra are then given by the Fourier 
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transforms of the corresponding current-current correlations ll45ll . 

/ + 00 
c?te*-*J||(g,t), 
'OO 

/+00 
dte''^'Mq,t). (4) 
-oo 

In simulations, Eqs. (51) are evaluated by using a discrete Fourier transform in a period of 1310 x 



III. THEORIES 



A. General formula 



A general expression for the dispersion relation, u;(q), in a magnetized 2D many-body system 



can be determined by solving the following equation for u; yi |5i 



|cu2l-M(q) -zC|| =0, 



where I is the 2D unit matrix, and i = ^/^. Note that the matrix 



29 



(5) 



UJLUr 



—IjJIjJr 







(6) 



accounts for the presence of an external magnetic field, with ojc = QB/mc being the gyrofre- 
quency of dust particles. Note that, although we have neglected damping in Eq. ([5]), it can easily 
be recovered by replacing cu'^ with lu^co + i'j). The matrix M(q) contains details of interactions 
between dust particles and its form depends on the specific theory used. 

For convenience, we introduce here a new parameter P = tOc/^pd- For typical laboratory dusty 
plasma parameters (uq = 1 x 10*^ cm^, ksTi = 0.1 eV, ksT^ = 3 eV, the radius of dust particles 
Td = O.lyum, m = 1.0 X 10"^^ g, Q = 100 e and n = 100 cm^), one has /3 = 0, 0.1, 0.5 and 1 for 
5 = 0, 1.5, 7.5, and 15 tesla, respectively. 



B. Random phase approximation (RPA) 



In the RPA, one obtains pOl] 



n 



MijpA(q) = — 0(A;)qq, 
m 



(7) 
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where (f){q) = 2ttQ a/ a/k^ + {qay is the 2D Fourier transform of the Yukawa potential ll30ll . The 
dispersion relation in the RPA can then be obtained by substituting Eq. ^ into Eq. ([5]) BOll as 



(8) 



where co'o(g) = uj'^^q^a? / a/ (ga)^ + k^. It is clear that the mode given in Eq. ([8]) is a hybrid 
cyclotron-plasmon mode in 2D resulting from a superposition of the cyclotron oscillations due to 
the magnetic field and a 2D dust-acoustic wave. It corresponds to the so-called magnetoplasmon 
mode in pure Coulomb systems , and is reduced to the usual 2D dust-acoustic mode when 
5 — i> [|8|]. More details of the derivation above can be found in Ref. [30]. 



C. Quasi-localized charge approximation (QLCA) 




FIG. 1: (Color onUne) (Color bar in arbitrary units) The wave spectra in both the j || q and j -L q directions 
with j3 = 0.5, K = 1 and F = 200. Dashed lines are the dispersion relations determined by the RPA and 
QLCA. 



In the QLCA, one has 

n 



MQLCA(q) = -0(A;)qq + D(q) 



m 



(9) 
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FIG. 2: (Color online) Profiles of the current-current correlation function in both the j || q (solid lines 
with filled circles) and j -L q (solid hnes with filled squares) directions at several fixed wavenumbers with 
/3 = 0.5, K = 1 and r = 200. 

The second term, Pfq) = iD ^uio)] (with fi and v = x, y), is the 2x2 dynamical matrix of the 



QLCA m, 



31 



46 



4811 . which takes into account the short-range inter-particle correlations 



that are missing in the RPA. By neglecting this term, we immediately recover the RPA result given 
in Eq. 



For the QLCA 113 ih . the wave dispersion relation has two branches: the high-frequency branch 
a;+(q) and the /ow-frequency branch t<;-(q), which are given by 
1 



± l\/[DLiq) + Driq) + uHq) + u^lf + AulDriq), 



(10) 



where Di^q) and Dxiq) are the projections of D(q) on the longitudinal and transverse directions, 
respectively. One recognizes that the high-frequency branch is actually the magnetoplasmon mode 
with corrections due to particle-particle correlations, while the low-frequency branch is identified 
as the magnetoshear mode, which is maintained by the strong particle-particle correlations. When 
B ^ 0, Eqs. (flOl) yield the longitudinal and transverse modes of the QLCA in unmagnetized 2D 



Yukawa liquids, respectively [|24l . 



28 



4l|] . It should be noted that the presence of the magnetic field 
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FIG. 3: (Color online) (Color bar in arbitrary units) The total wave spectra (sum of the j || q and j -L q 
directions) for /3 = 0.5, k = I and F = 200. Dashed lines are the dispersion relations determined by the 
RPA and the QLCA. 

does not alter the general scenario of the longitudinal and transverse modes at 5 = 0, although it 
strongly modifies the dispersion relation. Namely, one sees that both Diiq) and Driq) appear in 
both the high- and low-frequency branches as a result of the interplay between the particle-particle 
correlations and the external magnetic field. In addition, let us mention that the magnetic field has 
no effect on the static structure of a system in equilibrium. Therefore, expressions for D r an d Dt 



are exactly the same as those obtained without the magnetic field [|24l 



express 
271,128, 



41, 



D. Harmonic approximation (HA) 



In the HA, it is assumed that the charged dust particles form a perfect triangular lattice in 2D 
and are located at the points of that lattice while performing thermal oscillations around their 
equilibrium positions. In this case, one needs to ado pt a n another strategy to obtain the interaction 
matrix and the dispersion relation. One finds Is, 2\a11 

d^ch ,^ 



280 



M 



HA 



;q) = -E 



(11) 



9 



where again n and u = x,y, and the summation over i includes all points on the triangular lattice 



21 



27 



2m 



The resulting dispersion relations correspond to the so-called magnetophonon modes [|30, and 
are quite similar to those of the QLCA. There are two branches, given by 



± 



1 



[M,. + Myy + + 4(M2 - M,,M, 



(12) 

where M^j^ is determined by Eq. (fTTI) . The wave propagation depends on the polarization angle 
9, and the corresponding dispersion relations are periodic functions of 9 with a period of n/6 due 



to the hexagonal symmetry I13i|2l|,|27i|28|]. The angular dependence arises because the system is 



now anisotropic at short wavelengths 
1.5 
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K=i r=io 
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FIG. 4: (Color online) (Color bar in arbitrary units) The total wave spectra (sum of the j || q and j -L q 
directions) for /3 = 0.5, k = I and T = 10. Dashed lines are the dispersion relations determined by the 
RPA and the QLCA. 



IV. RESULTS AND DISCUSSIONS 

As a typical example, we show in Fig. [T] the simulation wave spectra in both the longitudinal 
(j II q) and transverse (j ± q) directions for {3 = 0.5, k = 1 and T = 200, along with the 
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FIG. 5: (Color online) (Color bar in arbitrary units) The total wave spectra (sum of the j || q and j -L q 
directions) for /3 = 0.5, k = I and F = 100. Dashed lines are the dispersion relations determined by the 
RPA and the QLCA. 

dispersion relations from the RPA and QLCA theories. The results clearly show that the high- 
and low-frequency branches co-exist in both the longitudinal and transverse directions. The high- 
frequency magnetoplasmon modes [|5D exhibit typical oscillations that originate from the short- 
range order in strongly coupled systems [|5D, whereas their long wavelength limits approach the 
cyclotron-plasmon oscillation frequency. However, the two branches of modes are not isotropic 
because the wave spectra appear to be different in the j || q and j -L q directions. This is shown 
particularly clearly in the low-frequency branches. For example, there is a gap at long wavelengths 
of the low-frequency branch in the j || q direction, and a gap around the first trough of oscillations 
in the j J- q direction. This indicates that coupling between the longitudinal and transverse modes 



strongly depends on both the wavenumber and frequency [12911 . A more detailed analysis of the 
frequency dependence of this coupling is shown in Fig.|2lfor several wavenumbers, where one can 
see that the amount of energy in the low- and high-frequency branches is redistributed between 
the two directions of wave propagation in a complementary manner. For example, when analyzing 
long wavelengths, one finds more energy in the low-frequency branch for the transverse direction 
than for the longitudinal direction, while there is more energy in the high-frequency branch for the 
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FIG. 6: (Color online) (Color bar in arbitrary units) The total wave spectra (sum of the j || q and j -L q 
directions) for /3 = 0.5, k = 1 and T = 1000. Dashed lines are the dispersion relations determined by the 
RPA and the QLCA. 

longitudinal direction than for the transverse direction. The roles of the transverse and longitudinal 
directions can be reversed at other wavelengths, as shown in Fig. [21 for qa = 3. We have observed 
in our simulations that such an effect occurs for any value of the coupling strength, and is more 
pronounced for higher /?s. We note that a similar effect was observed previously in the MD 



simulation of a perfect triangular lattice n2% 



In the following, we shall simply add the wave spectra obtained in both directions to achieve 
a more complete picture of both the high- and low-frequency branches. This is done in Fig. [3] 
using the data from Fig. [U where both branches now show similar oscillatory structures. One 
observes that the QLCA captures the dispersion of the high-frequency wave spectra very well in 
a reasonably large range of wavelengths, whereas the RPA displays noticeable discrepancy with 
our simulations except in the long wavelength limit. For the low-frequency branch, the agreement 
between the QLCA and simulations is less satisfactory, although the QLCA seems to give the 
correct trend in the oscillatory structures of the spectra. 

We first examine how collective dynamics depends on the state of the system. Figures |4l [51 and 
[6lshow the wave spectra for a high temperature liquid with F = 10, a cool liquid with F = 100, 
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FIG. 7: (Color online) (Color bar in arbitrary units) The total wave spectra (sum of the j || q and j -L q 
directions) for f3 = 0.5, k = I and T = 1000. Solid lines are the angle-dependent dispersion relations 
determined by the HA, with light-gray (yellow) for the high-frequency branch and dark-gray (red) for the 
low-frequency branch, respectively. 

and a crystalline state with T = 1000, respectively, where k = 1 and P = 0.5 are kept fixed 
in all three figures. One sees from Fig. |4]that, at F = 10, both the high- and low-frequency 
modes are heavily damped at short wavelengths, say qa > 1, and there are essentially no collec- 
tive modes beyond qa = 6. One would attribute this to Landau damping and viscous/coUisional 
damping. The QLCA dispersion relation again agrees reasonably well with our simulations for 
the high-frequency branch, in the region where the peak in the spectra is distinguishable, while the 
agreement between the QLCA and our simulations is poorer for the low-frequency branch. There 
may be several reasons for this. On one hand, since the low-frequency mode is always heavily 
damped in the simulations, the corresponding wave spectra spread out in the uj — q plane, and it 
is difficult to isolate their peak positions. On the other hand, it is well known that most of the 
damping effects, such as Landau damping and diffusional damping, are not included in the QLCA 



24 



m 



Am . This can be related to the inability of the QLCA to predict the long wavelength 
cutoff in the shear mode, which is clearly seen in the low-frequency branch in Fig. HI 

As expected, the damping effect becomes much weaker with increasing F or, equivalently, with 
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decreasing temperature. When T = 100, one observes from Fig.[5]that the wave spectra are largely 
extended in the short wavelength direction, and so is the region of agreement between the QLCA 
and the simulations for the high-frequency branch. However, the discrepancy between the QLCA 
and the simulations for the low-frequency branch is still significant. A substantial improvement is 
observed only in Fig. [6] at T = 1000, for which the long wavelength cutoff in the low-frequency 
branch disappears and oscillatory structures are quite distinguishable in both branches. A good 
agreement between the QLCA and simulations is found in both branches, particularly in the first 
Brillouin zone. For shorter wavelengths, some fine structures appear in the spectra due to the 
angular dependence of the wave dispersion, and the QLCA dispersion can only give a correct 
tendency of the oscillations. The task of resolving full angular dependence of the wave dispersion 
lies beyond the capability of the QLCA, but this issue can be tackled by the HA, as shown next. 




FIG. 8: (Color online) (Color bar in arbitrary units) The total wave spectra (sum of the j || q and j -L q 
directions) for (3 = 0.3, k = I and F = 1000. Dashed lines are the dispersion relations determined by 
the RPA and the QLCA. Solid lines are the angle-dependent dispersion relations determined by the HA, 
with light-gray (yellow) for the high-frequency branch and dark-gray (red) for the low-frequency branch, 
respectively. 

Figure |7] displays the same spectra as shown in Fig. [6l but comparisons are made with the 
dispersion results from the HA. There are two groups of dispersions from the HA, corresponding 
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to the high- and low-frequency modes, respectively, and in each group ten different curves that 
cover the period of 7r/6 evenly are shown. One sees that the structure of the wave spectra is well 
overlapped by the HA dispersion curves, indicating a good agreement between the HA and our 
simulation. By comparing the results of the HA and those of the QLCA in Fig.[6l one concludes 
that the QLCA gives a kind of angle- averaged dispersion. 




0.5 1.0 1.5 . 




FIG. 9: (Color online) (Color bar in arbitrary units) The total wave spectra (sum of the j || q and j -L q 
directions) for /3 = 1.0, k = I and F = 1000. Solid lines are the angle-dependent dispersion relations 
determined by the HA, with light-gray (yellow) for the high-frequency branch and dark-gray (red) for the 
low-frequency branch, respectively. 

Next, we analyze the dependence of wave spectra on the intensity of the external magnetic 
field. Figure [8] shows wave spectra for P = 0.3 (with other parameters being the same as those 
in Fig. |7]), together with the dispersions from both the QLCA and the HA. Comparing with Fig. 
Ill one sees from Fig. [8] that the oscillations in the wave spectra are greatly enhanced in this case, 
and fine structures due to the angular dependence in the wave propagation direction are enriched. 
As before, the agreement with the QLCA is clearly shown for both branches, roughly within the 
first Brillouin zone. On the other hand, one sees that the HA captures the fine structures in the 
wave spectra very well. The fact that the wave spectra in the high-frequency branch exhibit a gap 
in certain directions of propagation at the wavelengths 6 < qa < 7 can be rationalized as follows. 
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FIG. 10: (Color online) (Color bar in arbitrary units) The total wave spectra (sum of the j || q and j -L q 
directions) for /3 = 0.5, k = 2 and T = 1000. Dashed lines are the dispersion relations determined by the 
RPA and the QLCA. 

Since our simulation of the crystal lattice unavoidably has defects and contains distinct domains 
with different orientations, the measured wave spectra are expected to exhibit some preference 
to certain angles of propagation based on the dominant crystal orientation. One can roughly es- 
timate such an effect by comparing the spectral distributions with the HA dispersion relations, 
as illustrated in Fig. [8l On the other hand, both the oscillations and fine structures in the wave 
spectra are suppressed for stronger magnetic fields, as shown in Fig.[9]for /3 = 1.0, along with the 
dispersions from the HA. One sees in that figure that the high-frequency branch is raised above 
^/^pd = /3 = 1.0 and its width is greatly reduced, as is the width of the low -frequency branch, 
when the magnetic field increases. 

Finally, we examine the effect of screening on the wave spectra. Figures [lOl [HI and Fig. [6] 
compare the wave spectra for screening parameters k = 2, 3 and 1, respectively, with V = 1000 
and P = 0.5 kept fixed. The dispersion relations from the RPA and QLCA are also shown for 
comparison. It is seen that, with increasing k, the oscillatory structures in the spectra gradually 
flatten out, so that the dispersions become almost featureless when k = 3.0. At the same time, the 
width of the low-frequency branch decreases continuously with increasing k and, for k = 3.0, the 
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FIG. 11: (Color online) (Color bar in arbitrary units) The total wave spectra (sum of the j || q and j -L q 
directions) for /3 = 0.5, k = 3 and T = 1000. Dashed lines are the dispersion relations determined by the 
RPA and the QLCA. 

spectra become very blurry due to heavy damping. 



V. CONCLUSIONS 



In this paper, we have used the Brownian dynamics simulation to generate wave spectra of 2D 
strongly coupled dusty plasmas in an external magnetic field. We have examined how these spectra 
are affected by the coupling strength in such systems and by the intensity of the magnetic field. In 
particular, we have compared the dispersion relations resulting from the simulation wave spectra 
with those obtained from several analytical theories, including the random-phase approximation, 
quasi-localized charge approximation and the harmonic approximation. Two branches of collec- 
tive modes are generally observed in our simulation: the high-frequency magnetoplasmon and the 
low-frequency magnetoshear modes. The RPA can only reproduce the high-frequency branch, 
i.e., the magnetoplasmon mode, and it gives a good approximation only in the long wavelength 
limit, but remains robust for increasing coupling strengths, even up to the crystalline state. The 
QLCA can reproduce both of the above modes, and it shows a particularly good agreement with 
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the simulation in the high-frequency branch for wide ranges of coupling strengths, screening pa- 
rameters and magnetic field magnitudes. However, the agreement between the QLCA dispersion 
relations and the simulation spectra is generally poor in the low-frequency branch, particularly in 
a high-temperature liquid state, probably because Landau and diffusional damping are neglected 
in the QLCA. The situation becomes better in a crystalline state where a good agreement between 
the QLCA and the simulation is found in the low-frequency branch within the entire first Brillouin 
zone. However, the QLCA cannot capture fine structures in the wave spectra due to angular- 
dependence of the wave propagation in the crystalline state. In this case, it is found that the HA 
can give a very good description of the collective dispersions observed in the simulation. 

Our simulations show that the wave spectra exhibit profiles at fixed wavenumbers that have 
variable widths in the frequency domain. These widths generally correspond to various dissipa- 
tion mechanisms that give rise to damping of the collective dynamics in strongly coupled systems. 
We note that we have deliberately chosen a very weak coUisional damping on the neutral-gas 
molecules in order to reveal the effects of the magnetic field and other parameters. While com- 
parisons of the observed widths cannot be made with the above analytical theories, we have found 
some general trends in our simulations. The strongest effect on the frequency widths of wave 
spectra was found to come from the coupling strength, with the widths strongly increasing with 
increasing temperature of the system. This is similar to the observations in unmagnetized SDCPs 
112 811 . With regards to the intensity of the magnetic field, the spectra were found to be narrower at 
higher magnetic fields, whereas lower magnetic fields are accompanied by broad, fine structures 
in the spectra reflecting a strong anisotropy in the direction of propagation. These structures are 
generally very well captured by the families of dispersion relations coming from the HA. Finally, 
we have found that increasing the screening parameter gives rise to narrower and flatter wave 
spectra, with the low-frequency branch reduced to a heavily-damped mode akin to low-frequency 
noise. Clearly, more theoretical work is necessary to fully understand the dissipation mechanisms 
in strongly-coupled dusty plasmas in the presence of a magnetic field. 
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